1 



CD 

in 



Physical Effects of Infrared Quark Eigenmodes in LQCD * 
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A truncated determinant algorithm is used to study the physical effects of the quark eigenmodes associated with 
eigenvalues below 400 MeV. This initial study focuses on coarse lattices (with O(a^) improved gauge action), light 
internal quark masses and large physical volumes. Four bellwether full QCD processes are discussed: topological 
charge distributions, the eta prime propagator, string breaking as observed in the static energy and the rho decay 
into two pions. 
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1. Basics 

In the truncated determinant approach (TDA) 
to full QCD, the quark determinant, ^{A) = 
det{H) = det{j^{I^{A) — to)) is spht-up gauge in- 
variantly into an infrared part and an ultraviolet 
parti. 



Table 1 

Coarse lattices studied. 



V{A)=Vir{A)Vuv{A) 



(1) 



The ultraviolet part, 'Duv, can be accurately fit 
by a linear combination of a small number of Wil- 
son loops ||]. The infrared part Vifj{A) is de- 
fined as the product of the lowest N\ positive 
and negative eigenvalues of H, with |Ai| < Aqcd 
(typically, ~ 400-500 MeV). The eigenvalues Xa 
of H are gauge invariant and measure quark off- 
shellness. The cutoff (for the separation in Eq. 
^ is tuned to include as much as possible of the 
important low-energy chiral physics of the un- 
quenched theory while leaving the fluctuations 
of In ViR of order unity after each sweep updat- 
ing all links with the pure gauge action (assuring 
a tolerable acceptance rate for the accept/reject 
stage) [1). This procedure works well even for 
kappa values arbitrarily close to kappa critical. 

Initial studies using TDA focus on the qualita- 
tive physical effects of the inclusion of the infrared 
quark eigenmodes. For this purpose, coarse lat- 
tices with large physical volumes are appropriate. 



Label 


Volume; Fermions; Gauge Action 


PW6 


64; n^^t = 840 K = .2180; 




Naive /3 = 4.5 


RW6 


64; nc„t = 840 K = .1950; 




0(a2) /3 = 6.8 


RW8 


84; ricut = 1780 K = .1950; 




0{a^) 13 = 6.8 



*TaIk presented by E. Eichten 



2. Bellwethers 

The coarse lattices given in Table |l| are be- 
ing studied on PC clusters. The O(a^) improved 
gauge action, P = 6.8[1.0{plaq) - 0.08268(rect) - 
0.01240(para)], was adjusted in Ref. Q to have 
approximately the same lattice scale a = OAfm 
as the naive gauge action at /? = 4.5. The physi- 
cal lattice size is 2.4 (fm) for PW6 and RW6 and 
3.2 (fm) for RW8; = .2190 for PW6, and .1960 
for RW6 and RW8; and eigenvalue cutoff scale is 
560 MeV for PW6 and RW6, and 445 MeV for 
RW8. In addition, 10^ x 20 lattices at /3 = 5.7, 
csw = 1-57 and n^ut = 520 {Ecut = 460MeU) 
are being generated on ACPMAPS. 

As shown in Figure Q it took about 10,000 full 
steps for the PW6 lattice configurations to equi- 
librate (reflecting critical slowing down - a few 
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Figure 1. Reaching equilibrium for the PW6 lat- 
tices. 



hundred suffice on small volumes). 

Four bellwethers can be used to characterize 
the physical differences between quenched and 
full QCD. They are discussed (in order of increas- 
ing difficulty to observe in lattice calculations) in 
the following four subsections. 

2.1. Topolopy 

The topological charge, Qtop, can be ex- 
pressed in terms of the eigenvalues of the Wilson- 
Dirac operator. 

This sum is quickly saturated by the low eigen- 
values. 

In full QCD configurations with very small 
eigenvalues of H are suppressed by the quark de- 
terminant factor. In particular, non-zero topo- 
logical charges must be suppressed in the chiral 
limit [mq — > 0). Furthermore, the functional de- 
pendence of the topological charge distribution, 
Pq, on the light quark mass niq is predicted by 
the chiral analysis of Leutwyler and SmilgaQ. 

Pq ^ Iqixf - Iq+i{x)Iq-,{x) (3) 

where 

x = l/2F/>2 ^Fm, <V^7A>, (4) 

Iq are modified Bessel Functions of order Q and 
V is the total space-time volume. 
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Figure 2. Eta prime propagator for the RW6 lat- 
tices. The valence term (squares), hairpins term 
(triangles) and total propagator (diamonds) are 
shown separately. 



The quantitative agreement with the expected 
behaviour of Eq. ^ has already been reported 
for the TDA method §. General agreement is 
also observed on the coarse lattices of the present 
studies. 

2.2. Eta Prime Mass 

The relation between the axial U{\) anomaly 
and the rj' mass is well understood in full QCD. 
For two light quarks {Nf = 2), — rri^ + uiq 
where Wg = 2Nfx/ f'^ and the topological sus- 
ceptibility is Vx =< Qtop >quonchod. The full 
rj propagator is the sum of a connected (valence 
quark) term and a disconnected (hairpin) term. 
Thus, in the continuum, the momentum space full 
propagator can be written: 

{p' + ml + ml)-^= (p2 + m2)-i - (5) 
mg(p2 + m2)-i(p2 + rtil + ml)-^ 

These separate terms and their sum are shown in 
Figure || for the RW6 lattices. The cancellation 
between the valence and hairpin terms in the full 
propagator is evident. 

2.3. Static Energy 

To date no convincing evidence for string 
breaking in full QCD has been presented using 
calculations of the static energy alone. However 
string breaking has been seen using the TDA 
method in 2D QED|§ and by the HMC method 
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Figure 3. Comparison of the static energy be- 
tween the RW6 lattices (sohd squares) and un- 
quenched 6^ (open circles) (and 8^ (solid circles)) 
lattices a.t (3 — 4.5. The static energy was ex- 
tracted from the ratio of time slices T = 2/T = 1. 



in 3D QCD|^. Studies of the cc and bb systens, 
lead to the expectation that virtual pair effects 
(below heavy-light meson pair production thresh- 
old) will soften the linear rise in the static energy, 
while above threshold the potential will flatten 
out (i.e.) the string will break. The heavy-light 
meson mass is 0.81 ± .02 for the RW6 lattices. 

Figure |3| shows the static energy for 200 RW6 
lattices versus the same number of unquenched 
6* lattices at /3 = 4.5. There is evidence of 
screening from the virtual pairs but no hard ev- 
idence of string breaking is found. Seeing string 
breaking will require more statistics (to study 
{T — 3)/{T ~ 2) with small error bars) and also 
the study of the RW8 lattices. 

2.4. Vector Meson Resonances 

For the RW6 and RW8 lattices at k = .1950, 
the rho mass is 1.33 and the pion mass is 0.205 
(in lattice units); hence the ir/p mass ratio is 
close to the physical value. For example, the rho 
propagator for the RW6 lattice is shown in Fig- 
ure ^. However, since this is a P wave coupling, 
the physical volume of the lattice must be large 
enough that the decay is allowed with the first 
nonzero momentum, Pmin = This requires a 
10^ lattice (RWIO) or creating a rho with initial 
momentum Pmm- Neither of these alternatives 
have been studied as yet. 



Figure 4. Rho propagator for the RW6 lattices. 



3. Status 

The present status of full QCD bellwethers is 
as follows: 

• The behaviour of the topological charge dis- 
tribution as a function of light quark 
mass m^Q - Seen. 

• The eta prime mass - m^, = + mg - 
Seen. 

• The static energy - string breaking. - In 
progress but needs more statistics. 

• Vector meson resonances - p tttt. - Yet 
to be studied in detail. 

Results for these four bellwether processes on 
coarse lattices should be available within a few 
months. 
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